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Abstract. We propose a new hierarchical method which 
uses dynamical arguments to find and describe substruc- 
tures in galaxy clusters. This method (hereafter h-method 
Q_i or h-analysis) uses a hierarchical clustering analysis to 
determine the relationship between galaxies according to 



their relative binding energies. We have tested from N- 



body simulations, the two following features of the pro- 
,—l posed method: 

1. It extracts subgroups which are much more stable dur- 
00 ing the cluster evolution than those given by other 
(~^) techniques. 

2. There exists a reasonable similarity between the struc- 



tures found when only the coordinates (x,y,v z ) pro- 



vided by "observations" are considered, and those 
found by using the six phase-space coordinates 

Q H We have applied this method to two Abell clusters: 
' ABCG151 and ABCG2670. Our results imply that 
5^ ABCG151 is separated into two clusters, one of them 
is again divided into two subclusters. 

ABCG2670 has no subclustering. Our method allows 
however to extract the most bound galaxies in its dynam- 
ical core. 
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1. Introduction 

We address here the question of how to extract subclusters 
- or substructures - from a galaxy cluster? 

The existence of substructure in galaxy clusters was 
first noticed from optical data by Baier (1977) and, in 
an extensive study, by Geller and Beers (1982). Since the 
work by Forman & Jones (1981) it has also been detected 
from X-ray observations. Such a small-scale structure may 



provide us with very valuable information on the inter- 
nal cluster dynamics, and/or on cosmological scenarios of 
cluster formation. Several works (in the hierarchical for- 
mation scenario) show in fact that the existence of sub- 
structure in a galaxy cluster could be related to intrinsic 
parameters as its evolutionary status, the internal amount 
and distribution of dark matter, etc. (see, e.g., Serna, Al- 
imi & Schol 1994, and references therein). 

Moreover, the interpretation of subclusters as recent 
merger events might explain some of the cD galaxy prop- 
erties, as their orientation with respect to the environment 
(West 1994a), and the peculiar velocity distribution in rich 
clusters (Merritt 1985; Tremaine 1990). Still in the hier- 
archical formation framework, it has been suggested that 
the frequency of substructure in galaxy clusters may con- 
strain the density parameter ft and the spectrum of pri- 
mordial density fluctuations (Richstone, Loeb & Turner 
1992; Lacey & Cole 1993; Kauffmann & White 1993). 

The development of methods to detect and identify 
subclusters is therefore necessary to link observations with 
theoretical predictions and models. A huge number of pa- 
pers has been devoted to this subject and several substruc- 
ture detection methods have been proposed in the past. 
An exhaustive list of these methods can be found, for in- 
stance, in Bird (1995) or West (1994b). We will just men- 
tion here some techniques not referred to by these authors, 
as the use of a correlation function analysis proposed by 
Salvador-Sole et al (1993), or the wavelet analysis on X- 
ray images by Slezak et al. (1994) and Grebenev et al. 
(1995). 

We could roughly summarize most of the substructure 
detection methods as essentially based on searching for 
overprobable coincidences and correlations in the space 
of positions and/or velocities. Most of them just inform 
that substructure is probably present (or not) in a galaxy 
cluster, or find those regions in the sky clearly dominated 
by subclusters. 
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We consider that the problem of identifying substruc- 
tures in a galaxy cluster cannot be separated of questions 
like: how can substructures be defined from a dynami- 
cal point of view?, are they transient or quasi-permanent 
phenomena?, is their evolution nearly independent on the 
whole cluster?. 

We present a new method to identify subclusters in 
galaxy clusters which is based on coherent dynamical ar- 
guments. The algorithm chosen to construct this method 
is that of hierarchical clustering methods. It is based on 
the mathematical theory of cluster analysis (see, e.g., An- 
derberg 1973). 

2. Hierarchical Clustering Method 

The clustering analysis provides different mathematical 
methods to distribute a set of objects into a set of clusters. 
Among these methods, hierarchical clustering gives the 
clearest insight into the cluster structure and allows for 
a relatively easy identification of cluster members. This 
kind of methods can be formulated in terms of operational 
concepts (e.g., Schorr 1976) and it was first applied in the 
astronomy field to study how the galaxies in a sample 
are gathered in groups or clusters by Materne (1978), and 
later by Tully (1980, 1987) and Gourgoulhon et al. (1992). 

2.1. General features of hierarchical methods 

A hierarchical clustering method starts with a sample of N 
objects (galaxies). The hierarchical relationship between 
these objects can then be obtained from the following pro- 
cedure (see fig.l): 

1) We initially consider that each object constitutes a 
group by itself, that is, we start by considering a collection 
of N groups: ({d}, {G 2 }, {(?„}). 

2) We define an affinity parameter s (for instance: dis- 
tances, inverse of forces,...), and an affinity or linkage func- 
tion S(s) (see below). We construct an array Sij giving the 
value of 5 for each couple of groups (G«, Gj). 

3) We then extract the two groups presenting 
the highest affinity: min(5y). These two groups Gk, 
Gi are then "merged" and replaced by a single 
group Gm = {Gk,Gi}. So we are left with N — 1 
groups: ({G 1 },{G kl },...,{G n - 1 }). 

4) The merging procedure is repeated until only one 
cluster contains all the TV objects of the sample : G\ is 
left. 

5) The merging sequence can be easily visualized un- 
der the form of a hierarchical tree (h-tree) . Such a graph 
shows from the smallest groups (or pairs) to the greatest 
structure including several groups (see Fig. 2). 



Start with N objects 



Construct the affinity array 
Eij 



search k,l so that Ekl is 
minimum 



Replace groups Gk, Gl by a single 
group Gkl 




Fig. 1. Flow diagram corresponding to the hierarchical method 
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2.2. Constructing a particular hierarchical method 

The hierarchical method described above can be applied 
to several kinds of subjects. In order to construct in a 



0.05 
X 



Fig. 2. Example of a hierarchical tree 
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coherent way a particular grouping method we need to 
specify: 

1. What does a substructure mean?. 

2. Which is the affinity parameter more closely associated 
to our substructure definition?. 

3. Which linkage method gives the best results? 

It is important to remark that the two last points cannot 
be separated from the first one. As a matter of fact, all 
the physical information specifying the nature of the par- 
ticular problem under consideration is contained in the 
definition of the affinity parameter. Different choices of 
S will therefore lead to different kinds of substructures 
in a sample. In the past, for instance, to extract groups 
and clusters from a sample of galaxies, Materne (1978) 
used the Ward criterion, which gives the increase, from 
two groups to their merged group, of some coordinate dis- 
persion measure; Tully (1980) used a quantity which is 
related to the inverse of the force between two groups; 
and Gourgoulhon et al. (1992) used instead the inverse 
of the merged group density. However, these authors did 
not try to interpret the physical meaning of their resulting 
structures. It is then difficult to understand, from a dy- 
namical point of view, the differences between the above 
methods and to decide which is the best suited one. 

2.3. Substructure Definition 

We then start by defining what a substructure means. 
Only provided with such a definition, we will be able to 
find and to test the best method to identify such a class 
of structures. 

Observational experience just gives an intuitive idea 
of what a substructure is. An objective definition is how- 
ever difficult to find. Simple definitions like "overdense 
regions" or "groups of galaxies with a mean spatial sep- 
aration smaller than their local environment" have sev- 
eral inconveniences. For instance, such simple definitions 
would identify as a substructure a group of galaxies which, 
from a mere coincidence, are passing close together at a 
given instant. 

In order to translate our intuitive idea to a more appro- 
priate theoretical definition, we note that an isolated rich 
cluster has a membership which does not change through 
its dynamical evolution, except for some few escaping par- 
ticles by collisional processes as ejection or evaporation. 
A system whose membership is considerably reduced in 
a relatively short time interval cannot be considered as 
a bound cluster. Such a situation is however much more 
complex for substructures within a galaxy cluster because 
they are not isolated. In principle, their membership could 
be modified through the system evolution by interactions 
with other substructures or by changes in the gravitational 
potential of the overall cluster. Nevertheless, the member- 
ship of a substructure must be much more "stable" than 
that of any random subset of galaxies without a particu- 
lar dynamical signification. The most strongly bound sub- 
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structures will have longer lifetimes than a weakly bound 
galaxy subgroup where the external influence is much 
more important. We have then decided to look for galaxy 
groups in a cluster with longer membership lifetimes than 
any random group. If they exist, we will call them "sub- 
structures" . 

2.4- The affinity parameter 

As we have quoted above, since all the physical informa- 
tion is contained in the affinity parameter, a good choice 
of s is crucial to obtain the greatest similarity between 
the obtained structures and those corresponding to our 
substructure definition. Since we are trying to find gravi- 
tationally bound groups, definitions of s just based on the 
galaxy positions are not appropriate because, even if the 
force between two spatially close galaxies is rather strong, 
they could be not mutually bound if their relative velocity 
is high enough to overcome their gravitational attraction. 

A measure of the boundness of a gravitational system 
is its relative binding energy. 

where 

fiij = mi rrijl (m, : + mj), (2) 

Obviously, a negative E t j for two galaxies within a galaxy 
cluster does not necessarily imply that they are bound. 
Processes like close collisions or external gravitational in- 
fluences could signify that they are unbound. Neverthe- 
less, the more negative the relative binding energy is, the 
more strongly bound are two galaxies against the above 
quoted collisional or external processes. Such a feature 
is observed in several N-body simulations. For instance, 
Kandrup (1994) noted that "gravitational N-body simu- 
lations evidence that particles initially with high binding 
energy tend to end up with high binding energy, low with 
low, etc." 

We then adopt such a quantity as our choice for s. In 
section 3.1 we will test, from N-body simulations, that an 
affinity parameter defined in terms of the relative bind- 
ing energies is in fact appropriate to outline substructures 
defined by their longer membership lifetimes. 

2.5. Linkage method 

The cluster analysis theory provides three different linkage 
techniques to construct the affinity array Sij for a given 
choice of s: 

1. Single linkage: each group is characterized by the 
longest s value needed to connect any group member 
to some other member of that group. 

Sij = min(£ M „) (3) 

veGj 
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2. Complete linkage: each group is characterized by the 
longest s value needed to connect every group member 
to every other member. 

Sij = max(E Miy ) (4) 

3. Centroid clustering: each merged group is character- 
ized by the s value between its "parent" group means 
or center of mass properties. 



(5) 

where Eq.q. is the affinity (relative energy) between 
the mean of i € and that of j € Gj . 

A comparative study of the main properties of these 
different linkage methods has been extensively performed 
in the literature (see, e.g., Anderberg 1973). The single- 
linkage method is one of the very few clustering techniques 
which can outline nonellipsoidal clusters, and the resulting 
h-tree is invariant to any monotonic transformation in the 
affinity parameter. It also has interesting connections with 
certain aspects of graph theory, as the problem of finding 
the minimum spanning tree. Its main disadvantage is that 
it can fail to distinguish two groups poorly separated in 
the s space. Complete linkage is also invariant to mono- 
tonic transformations of the affinity measure. In contrast 
to the single-linkage method, it can only be interpreted in 
terms of the relationships within individual groups, while 
the differences between groups are not determined very 
reliably. Although it has some common aspects to a maxi- 
mally connected subgraph, it has no special interpretation 
in graph theory terms. Finally, the centroid method and 
its variants has the great disadvantage that reversals can 
occur, i. e., the Sij value corresponding to the groups to 
be merged may rise and fall from stage to stage of the hier- 
archical "merging" sequence. When these reversals occur, 
the resulting h-tree has no meaningful interpretation. 

In the next Section, we will also test from N-body sim- 
ulations each of these three linkage methods. We will show 
that structures within a galaxy cluster are much better 
described by the single linkage clustering given by Eq.(3). 
We then adopt such a linkage technique to fully specify the 
particular hierarchical clustering method that we propose 
here. 

3. Tests from N-Body Simulations 

In this section, we will show the efficiency of our method 
by testing: 

— The ability of a clustering method to extract those 
galaxy groups with the highest lifetimes. Different dy- 
namical arguments seem to indicate that such a pro- 
cedure is in fact appropriate to do it. However, it is 
necessary to test more rigorously this point not only 
by comparing with other possible choices of the affin- 
ity parameter s, but also by comparing the three main 
linkage techniques. 



— The ability of the method applied on observational 
data (position on the sky and line-of-sight velocity 
i.e. 3-D-data) to recover the structures obtained us- 
ing (6-D)-data. The efficiency of a clustering method 
must also be measured in terms of its applicability to 
real data, where some coordinates are unavailable. A 
method giving very different results when only the co- 
ordinates provided by observations are used, cannot be 
considered as useful in actual practice. 

The first point will be analyzed by performing N-body 
simulations of the self-gravitating evolution of a galaxy 
system, while the second one will require a series of static 
numerical experiences. 

In order to evaluate the stability of results when we 
compare, for a given system, two different evolutionary 
stages, or 6-D against 3-D configurations, we have im- 
plemented the following procedure. We start by applying 
the method under consideration to obtain the hierarchical 
trees (see Fig. 2) corresponding to each of the two con- 
figurations to be compared. For each of both trees, we 
construct an array defined so that the element is the 
number of 'mergers' in the hierarchical sequence involv- 
ing some group which contains simultaneously both i and 
j galaxies. The number of coincidences between the two 
trees (denoted by superindexes 1 and 2) is then 



i j^ti 

while that of differences is 



N, 



(6) 



(7) 



In order to quantify how different or similar are these two 
hierarchical trees, we then define the rate of instability or 
difference between them as 



S 



N, 



N c + N d 



(8) 



Obviously, the most stable method will lead to the small- 
est <S values. Two exactly equal trees imply 5 = 0, while 
those absolutely different lead to S = 1. 

3.1. Stability during the system evolution 

We have considered in our N-body simulations a rather 
non-stationary system without very high density con- 
trasts for substructures. Initial conditions were gener- 
ated by considering 400 equal-mass particles with random 
positions within an uniform sphere and radial velocities 
(spherical collapse model). The self-gravitating evolution 
of that system was then simulated by using the N-body 
code developed by Alimi & Scholl (1993) on Connection 
Machine. At t\ = t co i/5, some very slight condensations 
begin to be apparent and we have taken this configura- 
tion as our starting reference system. This configuration 
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will then be compared to the later one at ti = 2t co i/5. 
Obviously, since substructures are not isolated and the 
overall gravitational potential is changing, the group mem- 
berships will vary somewhat from ti to t^. We will now 
test the stability of structures defined by different meth- 
ods during the system evolution. 

3.1.1. The stability of linkage types 

Fig. 3a shows the resulting stability parameters S when 
each of the three linkage types (single, complete and cen- 
troid) are used together with the definition Eq.(l) for the 
affinity parameter. We see that the smallest evolutionary 
variations of substructures are obtained when a single link- 
age technique is used. Complete linkage leads instead to a 
S value greater by a factor of ~ 2 during the considered 
time interval, and centroid clustering leads to a even worse 
evolution stability. 

This result shows that the single-linkage technique cho- 
sen (by us) in our hierarchical clustering method is indeed 
the best suited to obtain stable substructures. This result 
can be interpreted by the numerical artifacts produced 
when complete or centroid clusterings are used. As we 
have quoted in Section 2.5, centroid and complete cluster- 
ing techniques are incapable of outlining non-ellipsoidal 
groups and they then introduce a spatial biasing in sub- 
structures. Moreover, if a centroid clustering is used, the 
value of Sij for two specific entities depends on the group 
sizes. This can result in a "snowball" effect where one of 
the first groups constructed by the hierarchical merging 
procedure starts to artificially accretate, due to its greater 
size, several single-component groups. 

3.1.2. The stability of different hierarchical methods. 

We have also tested the stability of our method (with a 
single linkage) as compared to other hierarchical methods 
proposed in the literature. Although some of these meth- 
ods were originally proposed to analyze large-scale cata- 
logs by using projected positions and, sometimes, also the 
line-of-sight velocities, their adaptation for substructures 
in a 6-D-simulated cluster is straightforward. 

Fig. 3b shows the stability parameter between ti and 
t<i obtained when different hierarchical methods are used. 
We see again that our binding method is much more sta- 
ble than other non-energy-based ones. In particular, the 
Ward criterion method proposed by Materne (1978) and 
also the density method by Gourgoulhon et al. (1992) 
are strongly unstable with the system evolution. Conse- 
quently, substructures defined by them do not correspond 
to bound subsystems in a cluster (except, presumably, if 
they are characterized by strong overdensities) . An inter- 
mediate stability is found however for the inverse-of-force 
method by Tully (1980). 

It is important to note that the good results obtained 
with our method are not simply due to the fact that it uses 



a larger amount of information (positions and velocities) 
than those just based on positions, like those involving 
forces or densities. The Ward criterion also uses velocities 
to measure the coordinate dispersion and, nevertheless, it 
leads to the worst results in terms of stability. We inter- 
pret this last result by the lack of dynamical coherence of 
such criterion. In the same way, one can expect that other 
dynamically non-coherent choices of the affinity parame- 
ter will also lead to unphysical results even if they use the 
same or much larger amounts of information. 

3.2. Stability against projections 

Each particle in a N-body simulation is characterized by 
its 6-D-coordinates {x, v}. However, astronomical obser- 
vations only provide 3 coordinates {a, S, v z }. We will now 
test the stability of different methods when only 3-D co- 
ordinates are used, that is, we will test if the " observed" 
structures are reasonably similar to the "intrinsic" (6-D) 
ones. 

The resulting stability parameters between t® and 
t\ D are shown in Fig. 3c (similar results are found for t\). 
We see that our method is the only one which does not im- 
ply considerably different results when projection biasing 
is present. The worst results are now obtained by Gour- 
goulhon's method, but Tully's and Materne's techniques 
also give very unstable results under projection effects. 
Concerning Tully's method, it must be noticed that, al- 
though it did not give very bad results when the evolution 
stability was examined, it implies however a rather poor 
similarity between observed and intrinsic substructures. 
This result clearly shows the efficiency of a hierarchical 
method based on relative energies as compared to other 
procedures. It is able to find relatively time-stable groups, 
and it also ensures a reasonable similarity between 3-D 
and 6-D structures. 

This behavior can be understood as follows: We quoted 
above that the structure of the h-tree obtained by our 
method is invariant to any monotonic transformation in 
the affinity parameter. This is roughly the case when our 
affinity parameter is transformed by projection effects. 
The relative binding energy is in fact the sum of two terms, 
the potential energy and the kinetic energy. When com- 
puted by using only the observed coordinates, both terms 
are separately smaller (or more negative) than the intrin- 
sic ones. Consequently, the "apparent" binding energy of 
any cluster galaxy will be smaller (more negative) than 
the 6-D one. However, since such effect is similar for all 
galaxies the relationships between them remain nearly un- 
changed. That is, two weakly bound galaxies close in the 
space of positions are usually also weakly bound under 
projection effects, as compared to the other ones, because 
their relative velocity holds rather high. 

It must be however noticed that, when real data com- 
ing from observations are used, the situation could appear 
more complicated because interlopers are often present. 
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For instance, let us consider the case of a galaxy with 
"reasonable" position and line-of-sight velocity , but with 
a very hight transverse velocity. Any method will (proba- 
bly) assign that galaxy to the "observed" cluster, while it 
does not actually belong to the 6-D-cluster. Such a situa- 
tion, although possible, can be assumed very rare in real 
clusters. 

The opposite case of a galaxy with a line-of-sight veloc- 
ity much higher (or much smaller) than the other galaxies 
in the field, will be easily interpreted by the h-analysis as 
an interloper, or as belonging to a different system, even if 
their apparent positions are close. It is precisely the case of 
the cluster ABCG151, we will discuss in the next section. 



4. Applications 

The method we propose, based on the relative energy and 
a single linkage, is efficient and stable in the various senses 
described above. We will now use this tool to analyze sub- 
structures in galaxy clusters. 



Fig. 3. Instability of hierarchical trees against dynamical evo- 
lution (a and b) and projection effects (c) for different linkage 
types (a) and different hierarchical methods (b and c) 



We transform angular distances between galaxies into 
distances in Mpc, by means of the cluster distance D = 
v^/H . H = 50km s -1 Mpc -1 is the Hubble constant, 
and ~vl is the mean redshift. Notice that when a clear par- 
tition is found after a first analysis, a further iterative 
study is performed by using the mean distances of sub- 
systems. We have used in the computation of the energy, 
a mass-to-light ratio for the overall cluster of the same 
order than that obtained by the virial theorem. After sev- 
eral computations we found that the resulting hierarchical 
trees for real clusters depend rather weakly on these pa- 
rameters, and they only indicate some differences if the 
mass-to-light ratio is over or underestimated by a factor 
of ~ 10 2 . 

We have processed ABCG151 and ABCG2670 with the 
h-method. 

The data for ABCG151 come from the ESO- 
Keyprogram on Nearby Abell Clusters of galaxies 
(ENACS Mazure et al 1991). The sample is described in 
den Hartog (1995). Notice that our data include all galax- 
ies in the direction of ABCG151, that is, we have not used 
redshift criteria to eliminate interloopers, foreground or 
background groups (see Fig. 4). This will allow us to show 
how such separation is performed by the h-method. 

Concerning ABCG2670, we have used the data given 
by Sharpies et al. (1988). In their paper, these authors 
provide us with a subsample of the galaxies belonging to 
the cluster (based on the redshift criteria) , and limited to 
magnitude bj = 19. We used only those galaxies located 
in a circle of radius r = 1200arcsec centered on the cD. 
This sample contains 48 galaxies. The velocity histogram 
is displayed in Fig. 6. 
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Fig. 4. Velocity histogram. Velocities are in km s 1 . The whole Fig. 5. Velocity histogram. Velocities are in km s . 
field of ABCG151, Gl and G2 (see text) ABCG151- Gl, G(l,l) and G(l,2) (see text) 



4.1. The cluster ABCG151 

We can clearly see on the h-tree shown in Fig. 7 that 
ABCG151 is constituted of two main clusters, Gl and 
G2. The method give explicitly those galaxies belonging 
to each of these clusters. Gl and G2 are superimposed 
on the sky, but they are separated in the "h-tree" (fig. 7) 
as well as in the velocity histogram (fig. 4). The mean ve- 
locity of Gl is V = 14896km s" 1 , while that of G2 is 
V = 29900km s _1 . If their velocities were only due to the 
Hubble expansion, the distance between the two clusters 
would be « 300Mpc 

Gl is subdivided into two subclusters G(l,l) and 
G(l,2). These two subclusters have nearly similar mean 
velocities as compared to the G2 velocity (see table 1). 
Clearly separated in the sky (Fig. 7), they perhaps actually 
form a pair of clusters. It is interesting to notice that these 
two clusters have velocity dispersions usual for clusters or 
large groups (see also Fig. 5). The analyses performed by 
the ENACS group, using fixed-gap and robust estimator 
for the velocity dispersion also lead to three clusters with 
the characteristics indicated in table 1 (bottom). 



Finally, G2 has no clear subclustering. The h-tree nev- 
ertheless shows a well bound pair of pairs in the deepest 
energy levels with some other neighbour galaxies forming 
a central dynamical core. A similar central structure will 
be found in the case of ABCG2670 discussed below (Sect. 
4.2). 

As we have explained above, we have also performed 
a h-analysis on Gl and G2 separately, but the resulting 
h-trees indicate minor changes only. 

4.2. The cluster ABCG2670: 

The first result we emphasize is that the h-tree does not 
reveal subclustering in ABCG2670 (Fig. 8). 

However more can be said: the deepest energy levels 
in the ABCG2670 h-tree show a well bound pair of pairs 
formed by galaxies (1, 19) and (6,8). It is interesting to 
remark that galaxy 1, the cD one, is located at the bottom 
of the energy well, and galaxy 19 is very close to it. The 
same behavior occurs for the pair (6,8) which is also lo- 
cated very close to the cluster center. These two pairs, to- 
gether with galaxies 44, 27, and 31 define a bound central 
group (denoted "b-c-g") in the deepest part of the energy 



8 



Arturo Serna and Daniel Gerbal: Dynamical Search for Substructures in Galaxy Clusters 




Fig. 7. A151 substructure analysis, a - the structural tree, b - the map 



levels. Again notice the similarity of central structures in 
ABCG2670 and ABCG151-G2. 

The velocity characteristics of ABCG27670 are dis- 
played in table 2. 

The result we found is very different from that given 
by C. Bird (1994) which is based on the analysis of the 
galaxy velocities and the hypothesis that subgroups have 
a Gaussian velocity distribution . However, this kind of 
hypotheses is questionable: we may in fact notice that the 
total luminosity of the cD galaxy is almost twice the to- 
tal luminosity of all the remaining "b-c-g" -galaxies. This 



feature suggests that these galaxies are rather a satellite 
system around gal.l instead of a virialized group with a 
Gaussian distribution of velocities. 

5. Comments on the hierarchical method 

Motivated by questions like: how can subclusters be de- 
fined from a dynamical point of view?; are they just tran- 
sient phenomena?, is their internal dynamics nearly in- 
dependent on the cluster where they are included?; we 
have developed a new method to identify substructures in 
galaxy clusters. 



Arturo Serna and Daniel Gerbal: Dynamical Search for Substructures in Galaxy Clusters 



9 




Fig. 8. A2670 substructure analysis, a - the structural tree, b - the map 



Since observational experience just provides us with an 
intuitive idea of what a substructure is, an objective an- 
swer to the first of the above questions is difficult to find. 
We have translated this intuitive idea to a theoretical defi- 
nition by looking for subsamples with longer lifetimes than 
any random subset of galaxies in a cluster. If they exist, 
we call them "substructures" . In order to implement such 
a definition, we have used a hierarchical clustering algo- 
rithm to determine the relationship between the cluster 
galaxies according to their relative binding energies. Ob- 



viously, a substructure definition different from what we 
propose could require a different working parameter. 

Since overdense regions in real clusters usually sat- 
isfy our above definition, the substructures defined by our 
method are not far from the intuitive idea of such kind of 
entities. Moreover, the inclusion of relative velocities in the 
kinetic term of binding energies (in addition to the spatial 
information contained in the potential term), results in a 
reasonable similarity between the observed and the intrin- 
sic substructures. In other words, the h-method is able to 
differentiate a cluster group from other groups or isolated 
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Fig. 6. Velocity histogram of ABCG2670. Velocities are in 
km s _1 . 





numbers 


mean 


median 


Vel. Disp. 


field 


105 


19612 


16310 




Gl 


72 


14896 


15515 


2230 


G2 


33 


29900 


30030 


846 


G(l,l) 


45 


15955 


15986 


748 


G(l,2) 


25 


12389 


12287 


390 




40 


29910 




857 




62 


15960 




667 




29 


12300 




395 



Table 1. Velocity characteristics of ABCG151. Top - from the 
hierarchical method - Field means: all the items in the field of 
the so-called "cluster". Gl, G2, G(l,l) and G(l,2). Bottom - 
ENACS analysis - Clusters are separated in the velocity space 
using a fixed gap. Velocities are in km s _1 . 





numbers 


mean 


median 


Vel. Disp. 


A2670 


48 


22744 


22646 


840 


cD 




23280 






"b-c-g" 


7 


22962 


22856 


571 



Table 2. Velocity characteristics of ABCG2670.b-g for the 
central bounded group, see text. Velocities are in km s _1 . 



galaxies just superimposed on the same region of the sky, 
to include galaxies dynamically bound to a group but with 
orbits or positions far from the region clearly covered by 
that group, and to identify subclusters not characterized 
by a remarkable overdensity in a cluster. All these fea- 
tures, as well as the evolution stability of structures found 
by our method, have been tested by means of numerical 
experiments and N-body simulations. 

Concerning the transitory character of substructures, 
the results obtained in Sec. 3.1 from N-body simulations 
show that their lifetime strongly depends on how they are 
defined and on the algorithm used to find them. Other 
previous methods as those referred to in Section 3.1, define 
substructures which could just be transient entities with 



rather short lifetimes. Our results show however that it 
is possible to define subclusters in such a way that they 
are much more stable throughout the cluster evolution. 
Groups located in the h-tree at the deepest energy levels 
will have a priori the longest lifetimes. Notice by the way 
that, since our method is able to extract substructures 
subsisting during long times, it appears specially useful for 
those works looking for substructures as a discriminant 
tool between various formation scenarios. 

Substructures are often considered in the literature as 
systems isolated from the whole cluster. We would like 
however to emphasize some aspects concerning this point: 

From a rigorous point of view, they can never be con- 
sidered as isolated because they are embedded in the clus- 
ter. Their dynamics are then influenced by other substruc- 
tures and also by the overall cluster. Nevertheless, by look- 
ing at the h-trees obtained for clusters like ABCG151 and 
AGCG2670, we see that some subgroups are characterized 
by a strong increase in their binding energy. This is for in- 
stance the case of galaxies (1,19,6,8) in the h-tree shown 
in Fig. 8. A priori, it can be expected that such a class of 
subsystems tend to be much more isolated, or indepen- 
dent, than other groups characterized by weaker increases 
of their energy levels. Our h-tree technique then provides 
an idea of how good is the approximation of "isolated" for 
a given galaxy group. 

Much work is still needed to improve the h-method. 
For instance, by including an estimation of confidence 
levels of the obtained substructures; by taking into ac- 
count the presence of a background non-luminous compo- 
nent (X-ray emitting gas and dark matter) dominating the 
gravitational potential. Although reasonable results could 
be expected if light traces the mass, a rigorous study of 
this problem is needed, specially in the opposite case with 
a considerable luminous-dark matter segregation. These 
improvements are in progress. 

Remarks: Information on the h-program can be pro- 
vided by e-mail at: serna@gin.obspm.fr 
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